Evidence for Superfluidity in a Resonantly Interacting Fermi Gas 
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We observe collective oscillations of a trapped, degenerate Fermi gas of 6 Li atoms at a magnetic 
field just above a Feshbach resonance, where the two-body physics does not support a bound state. 
The gas exhibits a radial breathing mode at a frequency of 2837(05) Hz, in excellent agreement 
with the frequency of vr = ^/l0v x v y /?> — 2830(20) Hz predicted for a hydrodynamic Fermi gas with 
unitarity limited interactions. The measured damping times and frequencies are inconsistent with 
predictions for both the collisionless mean field regime and for collisional hydrodynamics. These 
observations provide the first evidence for superfluid hydrodynamics in a resonantly interacting 
Fermi gas. 
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Strongly-interacting two-component Fermi gases pro- 
vide a unique testing ground for the theories of exotic 
systems in nature, ranging from super-high tempera- 
ture superconductors to neutron stars and nuclear mat- 
ter. The feature which all of these systems have in 
common is a strong interaction between pairs of spin- 
up and spin-down particles. In atomic Fermi gases, 
tunable, strong interactions are produced using a Fes- 
hbach resonance [H 0- Q Near the resonance, the 
zero-energy s-wave scattering length a exceeds the in- 
terparticle spacing, and the interparticle interactions are 
unitarity- limited and universal 0, 0-0]. In this re- 
gion, high temperature Cooper pairing has been pre- 
dicted [aSHEi 

In this Letter, we present measurements of the fre- 
quencies and damping times for the radial hydrodynamic 
breathing mode of a trapped, highly degenerate gas of 
6 Li atoms just above a Feshbach resonance These 
observations provide the first evidence of superfluid hy- 
drodynamics in a resonantly interacting two-component 
Fermi gas. 

Hydrodynamic behavior in a collisionless quantum gas 
at very low temperature is known to be a hallmark of 
superfluidity. Previously, we observed hydrodynamic, 
anisotropic expansion of a strongly interacting, ultracold, 
two-component Fermi gas [l^- However, an initially col- 
lisionless gas could have become collisionally hydrody- 
namic as the Fermi surface significantly deformed during 
the expansion |13|. Thus, the observations suggested su- 
perfluid hydrodynamics, but were not conclusive. 

Superfluidity has been observed in Bose-Einstein con- 
densates (BECs) of molecular dimers, which have been 
produced from a two-component strongly interacting 
Fermi gas. The first experiments produced the dimcr 
BECs at a magnetic field below the Feshbach reso- 
nance, where the atoms have a nonzero binding en- 
ergy 0, 0, 0, 0. Although the two-body physics 
does not support a bound state at magnetic fields above 
the Feshbach resonance, for fermionic atoms, the many- 
body physics does. Observations of BECs originating 



from such dimers are consistent with the existence of 
preformed pairs |l8l Il9| . These experiments indirectly 
explore the microscopic structure, while our experiments 
are complementary in that they directly measure the 
macroscopic dynamics. 

One method for distinguishing between a BEC and a 
superfluid Fermi gas is to examine their collective hydro- 
dynamic modes at low temperature where the trapped 
gas is collisionless [2(| ■ For a weakly repulsive BEC con- 
tained in a nearly cylindrically symmetric trap, the radial 
breathing mode occurs at a frequency of 

v B = 2^i/ x v y , (1) 

where v% are the harmonic oscillation frequencies in Hz 
of a noninteracting gas in the i th direction of the trap. In 
contrast to a weakly repulsive BEC, a superfluid Fermi 
gas in a cigar-shaped trap is predicted to have a radial 
breathing mode at the hydrodynamic frequency 



VH = V ~?~ VxVy ' ( 2 ) 

This result is obtained in the unitarity limit, where the 
shift from the interparticle interactions vanishes for a hy- 
drodynamic gas |2lj,|22|. In general, hydrodynamics with 
the frequency uh can arise from superfluidity or from col- 
lisions in a normal fluid. However, at low temperature, 
Pauli blocking is expected to suppress the collision rate 
in a Fermi gas, increasing the damping of the collisionally 
hydrodynamic modes as the temperature is lowered. 

In our experiments with a trapped Fermi gas, Pauli- 
blocking of collisions is expected to be effective at the 
lowest temperatures achieved. Nevertheless, a weakly 
damped radial mode at precisely the hydrodynamic fre- 
quency of vu is observed, and the damping rate decreases 
strongly as the temperature is lowered below ~ 30% of 
the Fermi temperature for a noninteracting gas. These 
observations provide the first evidence for superfluid hy- 
drodynamics in a resonantly interacting Fermi gas. 

We measure the frequencies and damping times of the 
radial breathing mode of the 6 Li Fermi gas near the Fes- 
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hbach resonance at 822 (3) G |l9(. The frequencies mea- 
sured at several temperatures and magnetic fields are 
compared to the predictions (Eqs.nand|2Jl based on the 
trap oscillation frequencies measured by parametric res- 
onance. A cross check of the measurement method is 
provided by measuring the breathing mode of the nonin- 
teracting gas. 

We begin by preparing a degenerate, 50-50 mixture of 
the lowest spin-up and spin-down states of 6 Li atoms in 
an ultrastable CO2 laser trap as described previously [lif . 
Forced evaporation at a chosen magnetic field in the 
range 770-910 G produces a highly degenerate, unitarity- 
limited sample. The trap depth is lowered by a factor of 
~ 580 over 4 s, then recompressed to 4.6% of the full trap 
depth in 1 s and held for 1 s to assure equilibrium. This 
produces a degenerate sample of N /2 ~ 1.5 x 10 5 atoms 
per spin-state at a temperature of ~ 0.15 Tp. The corre- 
sponding Fermi temperature for a noninteracting gas is 
then Tp = fiui^N) 1 ^ = 2.5 (u> = 2ir x {v x v v v x )V 3 ), 
small compared to the final trap depth of 35 fiK. 

From the images of the released cloud, the number of 
atoms is determined by integrating the column density. 
We have measured the resolution of the entire imaging 
system to be 5.5 /im. For low temperatures, T/Tp < 0.4, 
the ratio T /Tp is determined by fitting a Thomas-Fermi 
profile for a noninteracting Fermi gas to the transverse 
(x) distribution obtained by integrating the column den- 
sity in the axial (z) direction. Very good fits are obtained 
for the noninteracting gas at 526 G. At 910 G, this proce- 
dure yields temperatures as low as T/Tp — 0.06 and ex- 
cellent fits. However, at fields closer to resonance, slightly 
higher temperatures are obtained, and the shape may 
not be precisely Thomas-Fermi due to many-body effects. 
For measurements near resonance at higher temperatures 
T/Tp = 0.5—1.2, where the gas is not perfectly hydrody- 
namic, temperature estimation is less precise. Here, we 
simply assume hydrodynamic expansion of a Maxwell- 
Boltzmann spatial distribution for a noninteracting gas. 
The measured temperatures therefore indicate the trend 
but not necessarily the absolute temperature. We find 
consistency between the measured number of trapped 
atoms, the temperature, and the initial cloud size ob- 
tained by hydrodynamic scaling 12]. For the strongly 
interacting gas, we include a reduction of the cloud ra- 
dius arising from the mean field ||. By correcting se- 
lected images for our estimated saturation I / I sa t = 0.2, 
we estimate that the true temperatures are lower by 0.03 
Tp and the true atom numbers are increased by a factor 
near 1.15 compared to the values given in Table H 

Trap oscillation frequencies at 4.6% of the full well 
depth are measured by parametric resonance in a weakly- 
interacting sample. The gas is cooled by forced evapora- 
tion over 25 seconds to temperatures of 0.3 Tp at a field 
of 300 G and the trap depth is then modulated by 0.5% 
for 1 s. During this period, the low collision rate pro- 
duces little damping, but permits the gas to thermalize. 



After modulation, imaging at 526 G is used to measure 
the release energy versus drive frequency. Well-resolved 
resonances are obtained at 2io x = 2ir x 3200(20)Hz and 
2u> y = 2tt x 3000(20) Hz. Because of the low frequency, 
the axial resonance is measured at full trap depth. We 
obtain 2lo z = 2ttx 600(20) Hz, yielding 2lo z = 2vrx 140(5) 
Hz at 4.6% of full trap depth, including a quadratically 
combined magnetic field curvature contribution of 21 Hz. 
From these measurements, v±_ = ^v x v y = 1550(20) Hz. 

To excite the transverse breathing mode, the trap is 
turned off abruptly (< 1 /is) and turned back on after a 
delay of to = 50 /is. To show that our excitation is a weak 
perturbation, we estimate the energy increase, which 
arises principally from the change in potential energy in 
the transverse directions. Assuming approximately bal- 
listic expansion, AE± = E^(cj^t ) 2 /2 = 0.1 Ej_. For 
initial temperatures of 0.1-0.15 Tp, the corresponding- 
temperature change is AT/T F ~ 0.05 when the gas ther- 
malizes, consistent with our measurements. To measure 
the frequency and damping time, the breathing mode is 
excited and the sample is held for a variable time thold, 
after which the trap is extinguished suddenly, releasing 
the gas. The cloud is allowed to expand for 1 ms and then 
imaged using a camera beam probe pulse of 5 /is dura- 
tion arranged to produce a two-level system as described 
previously |12(. 

To study a noninteracting sample, breathing modes 
are excited at 526 G, where the scattering length is 
nearly zero |U 0]> after cooling at 300 G as de- 
scribed above. Fig. ^ shows y/ (x 2 ) for the expanded 
gas ((x) = 0), plotted as a function of the hold time 
in the trap, thold- Fitting an exponentially damped sinu- 
soid x rms + Aexp(— i/Vdamp) sin(27rz/t -I- tp) to the data, 
we obtain the frequency v — 3212(30) Hz and the damp- 
ing time Tdamp = 2.04 ms, where the errors are from the 
fit. The damping time is consistent with a small anhar- 
monicity of the gaussian profile of the trap potential. 




thold( m s) 

FIG. 1: Excitation of the breathing mode in a noninteracting 
Fermi gas of 6 Li at 526 G. Error bars denote one-cr uncer- 
tainty. 

The strongly interacting gas exhibits longer damping 
times at low temperature. Fig. |2Ia-c) shows the results 
at 870 G as the temperature is lowered from T/Tp — 0.50 
to 0.17. The damping time increases from 1.4(0.1) ms to 
3.85 (0.4) ms. At the lowest temperature, the measured 
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FIG. 2: Excitation of the breathing mode in a strongly in- 
teracting Fermi gas of 6 Li at 870 G. a) T/T F = 0.50, b) 
T/T F = 0.33 c) T/T F = 0.17. Error bars denote one-cr un- 
certainty. 



oscillation frequency is 2837(05) Hz. 

Results for different temperatures and magnetic fields 
are given in Tabled 

The measured oscillation frequencies for the transverse 
breathing mode can be compared to those expected for a 
noninteracting gas, for a weakly repulsive Bose gas, and 
for a hydrodynamic Fermi gas. For a weakly repulsive 
Bose gas, using eqs. ^ and |21 we predict vb = 2z/j_ = 
3100 Hz and for a hydrodynamic gas, vh = y/lQ/3vj_ = 
2830(20) Hz. 

The measured oscillation frequency of 3212(30) Hz for 
the breathing mode of the noninteracting gas is in excel- 
lent agreement with the frequency 3200(20) Hz measured 
by the parametric resonance method for the x-direction 
which is imaged in the experiments. Hence, the paramet- 
ric resonance method and the breathing mode excitation 
method yield identical frequencies to better than 1%. 

For the strongly interacting gas at 870 G and T/T F = 
0.17, the measured radial breathing mode frequency 
2837(05) Hz, is in excellent agreement with the predic- 
tion vh = 2830(20) for a hydrodynamic Fermi gas, and 
it differs significantly from that of the noninteracting gas 
and the weakly interacting Bose gas. Similar results are 
obtained at 860 G and 880 G and 910 G at slightly lower 
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2836(16) 


3.95(1.5) 
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2798(15) 
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TABLE I: Breathing mode frequencies v and damping times 
Tdamp- B is the applied magnetic field, T/Tf is the initial 
temperature in units of the Fermi temperature and N is the 
total number of atoms, uncorrected for saturation. a; rms is the 
time-averaged root-mean-square size of the oscillating cloud. 
^From the tails of a bimodal distribution. * Shown in the 
figures. L For t = 25 fis. 2 At 18.8(0.9)% trap depth, to = 
25 fjs and 0.8 ms expansion time. Error estimates are from 
the fit only. 



temperatures. For the axial direction, the measured am- 
plitude of the oscillation is consistent with zero. 

On the molecular side, just below resonance at 815 
G at T/T F = 0.14, we obtain v = 2931(19), which is 
just above the vjj, consistent with predictions that the 
response near resonance is fermionic • At a much 

lower field of 770 G, we find v = 3000(150) closer to the 
predicted Bose frequency of 3100 (20) Hz. However, the 
data is not of as high quality as that shown in Fig. [2J 

Over the range of magnetic fields studied, our mea- 
sured oscillation frequencies v (B) at the lowest tempera- 
tures show the same ma gne tic field dependence as those 
of Rudi Grimm's group 25] . However, our data show a 
much smaller frequency shift with respect to the expected 
hydrodynamic frequency vr- 

We estimate the frequency shifts Av = v(me&s) — 
^(actual) arising from anharmonicity in the trapping po- 
tential |26J. For the hydrodynamic frequency, Ai/jj = 
-(32/25) y/l0/3i> ± Mu)l_x? m J(b 2 H U), where U is the 
trap depth, M is the 6 Li mass, and uj± = 2irv^. Here, 
bu is the hydrodynamic expansion factor, 11.3 after 1 
ms jiij ■ The shift in the geometric mean of the trans- 
verse frequencies Av±_ = —( t o/5)v±Mu> 2 1 _x^ rlls /(b B U) 1 
where b B — 10.3 is the ballistic expansion factor at 
1 ms. Using Table these results yield a net Av = 
Av H -y/T0/3Av ± = +24 Hz at 870 G and T/T F = 0.17. 
For the three higher temperatures at 870 G, we find 
Av ~ -35 Hz at T/T F = 0.33 and 0.5, and -16 Hz at 
T/T F = 1.15, consistent with the measured ~ —60 Hz 
shift below the lowest temperature data. 

We have also investigated the effect of decreasing the 
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FIG. 3: Product of the damping time and the measured 
breathing mode frequency versus temperature. The open cir- 
cle shows the result when the trap depth is increased by a 
factor of 4.1. 



oscillation amplitude to less than 10% by reducing to to 
25/iS at 870 G and T/T F = 0.15. For small amplitudes, 
the aspect ratio of the cloud changes very little. Hence, 
we expect that the deformation of the Fermi surface is 
very small, so that collisional behavior is not induced. 
We obtain v = 2838(16) Hz, in precise agreement with 
the results for to — 50 fis. The damping time is somewhat 
increased to 6.00(1.4) ms, presumably due to the smaller 
energy input, rather than reduced anharmonicity (since 
the frequency is unchanged). 

The measured damping time shows a rapid increase 
with decreasing temperature, Fig. The open circle 
shows the result at a trap depth increased by a factor 
of 4.1. The measured hydrodynamic frequency scales 
as a/4-1 within 3% and the product v meas Tdamp is con- 
sistent with those at lower trap depth. This is consis- 
tent with a damping rate which scales linearly with trap 
frequency, as expected for unitarity-limited interactions, 
where the rate scales with the Fermi energy. Note that 
anharmonicity cannot make a major contribution to the 
temperature dependence of our damping rates: The an- 
harmonic contribution to 1/ Tdamp would be proportional 
to the shift and therefore independent of trap depth for 
fixed T/Tp and N. Then, vTdamp °c and the vertical 
position of the open circle would increase by a factor of 2. 
Further, the anharmonic contribution is identical for the 
three highest temperature points where the mean cloud 
sizes are nearly the same. 

We have attempted to model the data at 870 G without 
invoking superfluidity. A first scenario is that the gas 
is nearly collisionless at the lowest temperature, and the 
long damping time and measured frequency are the result 
of collisionless mean field evolution. A second scenario is 
collisional hydrodynamics. 

The collisionless mean field scenario requires a large 
negative mean field shift to explain the difference be- 
tween the frequencies of 3212(30) Hz and 2837(20) Hz 
measured for the noninteracting and strongly interacting 



samples, respectively. However, for a unitarity limited 
interparticle interaction with a negative j3 = —0.55 [28|. 
a Vlasov equation model [2^ yields a +90 Hz shift rela- 
tive to 3200 Hz, while we observe a -400 Hz shift. Also, 
for our trap, the same model shows that the coupling 
of the collisionless transverse modes by the interaction 
would produce a noticeable beat at 370 Hz with an am- 
plitude minimum at 1 ms, which is not observed. Hence, 
the data are inconsistent with the collisionless scenario. 

To investigate the second scenario, we considered a col- 
lisional h ydr odynamic model, including two-body Pauli 
blocking [2fJ, |30|. A small negative shift at the higher 
temperatures might arise from the mean field in the hy- 
drodynamic limit [2lL E2T | , or from the anharmonic shift 
described above. Ne glec ting these shifts, a relaxation 
approximation model |3jJ can be used to determine both 
the breathing mode frequency and the damping time in 
terms of the measured trap oscillation frequencies, for an 
arbitrary momentum relaxation rate. We find that a very 
large momentum relaxation rate is needed to fit the 4 ms 
damping time of the T/Tp — 0.17 data in a collision- 
ally hydrodynamic regime. Then, the predicted damping 
time is large over a broad temperature range, inconsistent 
with the observed rapid decrease in damping time with 
temperature. Lowering the maximum relaxation rate, 
we can fit the damping times at the two highest temper- 
atures. In this case, however, obtaining a 4 ms damping 
time requires a temperature below T/Tp = 0.1, i.e., a 
nearly collisionless regime, inconsistent with observations 
as described above. 

In conclusion, at our lowest temperatures, we observe a 
breathing mode at precisely the hydrodynamic frequency 
as well as highly anisotropic hydrodynamic expansion, as 
in our previous experiments 12, 29]. The damping time 
increases rapidly as the temperature is lowered [32j , con- 
sistent with a transition from collisional to superfluid hy- 
drodynamics at a temperature between 0.2 and 0.3 Tp. 
On the basis of the above arguments, we believe the data 
are not consistent with either collisionless mean field evo- 
lution or collisional hydrodynamics. It is therefore diffi- 
cult to see how the observations can be explained without 
invoking superfluidity. 

Recent theory describes the BCS-BEC crossover 
regime in terms of very large fermionic pairs, compa- 
rable in size to the interparticle spacing [33I I34I 
Falco and Stoof 34] predict BEC-like or BCS-like behav- 
ior above the Feshbach resonance, depending on whether 
e b = fi 2 1 '(a 2 m at0 m) is < 2k B T F , or > 2k B T F j^j. Near 
resonance, where e& << 2kgTp, the majority of fermionic 
pairs (either Bose molecules or Cooper pairs) are very 
large. Hence, one expects that the response of the sys- 
tem to compression is fermionic, and scales with density 
as ri 2 / 3 1^3]. This idea is consistent with our measure- 
ments and with predictions [2lL Eij ] . 
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